AD 624703. 


Operations Research and Economics Division 


SENSITIVITY ANALYSIS OF 
CIVIL-DEFENSE ^SYSTEMS AND COMPONENTS 

FINAL REPORT VOLUME M 
R-OU-157 


A GENERALIZED SENSITIVITY ANALYSIS OF CD 
SYSTEMS 



Prepared for 


Office of Civil Defense 
United States Department of the Army 
under 


Office of Civil 

dtoSCQMoct) 


Defense Contract No. 
OCD Subtask 4113E 
RTI Project OU-157 


OCD-PS-64-56 







THE RESEARCH TRIANGLE INSTITUTE 


Operations Research and Economics Division 


PCD Review Notice 


This report has been reviewed in the 
Office of Civil Defense and approved 
for publication. Approval does not 
signify that the contents necessarily 
reflect the views and policies of the 
Office of Civil Defense. 

Distribution of this document is 
unlimited. 


RESEARCH TRIANGLE INSTITUTE 
Durham, North Carolina 

FINAL REPORT VOLUME III 
R-OU-157 


A Generalized Sensitivity Analysis of CD Systems 
by 


John H. Neblett 


and 


Kenneth E. Willis 
October 1, 1965,. 

Prepared for 



Office Of Civil Defense 
United States Department of the Army 
under 

Office of Civil Defense Contract No. OCD-PS-64-56 
OCD Subtask 4113E 
RTI Project OU-157 


VD€5^403 







FINAL REPORT VOLUME III 


R-OU-157 


A Generalized Sensitivity Analysis of CD Systems 


Prepared for 

Office of Civil Defense 
Unites States Department of the Army 
under 

Office of Civil Defense Contract No. OCD-PS-64-56 
OCD Subtask 4113E 
RTI Project OU-157 

by 

John H. Neblett 
and 

Kenneth E. Willis 


RESEARCH TRIANGLE INSTITUTE 
Operations Research and Economics Division 
Post Office Box ''90 
Durham, North Carolina 



Approved by: 



Edgar A 


j2 




. Parsons, Director 


R. S. Titchen, Deputy Director 


1 October 1965 






FOREWORD 


This is Volume III of three separately bound volumes in which are reported th 
research completed under the general terms of the Office of Civil Defense Subtask 
Number 4113E, "Sensitivity Analysis of Civil Defense Systems and Components." 

The research of the authors was very ably supported by Mr. Herbert Hill, and 
Miss Mary B. Woodside. Mr. Hill assisted in the development of the models and 
Miss Woodside performed the calculations reported in several of the Appendixes. 


iii 








ABSTRACT 


This sensitivity analysis employs mathematical models which estimate the total 
dose. Maximum ERD, and probability of casualty or fatality for an individual exposed 
to a particular radiological environment. The objective of the analysis is to 
determine the relative importance of the parameters normally employed in evaluating 
Civil Defense systems; i.e., to determine those parameters which, through large 
variance or inaccurate estimates, will contribute most to erroneous evaluations of 
CD systems. The sensitivity of the dose (or probability of casualty) to variations 
in the input parameters defining the radiological environment is examined. The 
total Cone model is analytical and the ERD model is computerized. The parameters 
examined are reference intensity, time of arrival, time outside in fallout, radia¬ 
tion decay exponent, ERD recovery fraction, ERD recovery rate, duration of fallout 
buildup, and protection factor. Sensitivity indices are calculated for each parame¬ 
ter. The sensitivity index is defined as the fractional change in dose (or proba¬ 
bility of casualty or fatality) divided by the corresponding fractional change in 
the input variable. It is concluded that dose and casualty computations are quite 
sensitive to errors in the field decay exponent, and they remain sensitive over 
the examined range. Sensitivity to variations in fallout reference intensity 
and protection factor are high over the whole range of parameter values. Sensitivity 
of time of arrival of fallout can be quite high for the lower values of the parame¬ 
ter. Sensitivity of dose and casualty computations to the remaining parameters is 
low in most cases of interest. Expansion of the sensitivity analysis to include 
parameters other than fallout, which define the total casualties from a given attack 
on the United States, is necessary before further conclusions concerning a national 
vulnerability analysis can be drawn. 
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A Generalized Sensitivity Analysis of CD Systems 


I. INTRODUCTION 

A. Purpose 

An analysis of the effectiveness of the civil defense shelter system must take 
into account many variables. These variables include some that determine the fall¬ 
out environment (fallout reference intensity, its time of arrival, its rate of build¬ 
up, and the radioactivity decay rate); some that determine the protection against 
radioactivity (the gamma ray shielding characteristics of available shelter, and the 
time required for travel through fallout before shelter is reached); and some that 
influence physiological condition (recovery rate, and the fraction of radiation 
damage that can be repaired). The research reported in this volume is directed at 
evaluating the importance of uncertainty or inaccuracy in each of these variables 
in estimating effectiveness of a shelter system in a fallout environment. Effective¬ 
ness is measured in terms of Equivalent Residual Dose (ERD), probability of fatality, 
and time required in primary shelter for the ERD to return to a specified level. 

The results from such an analysis can be applied in two ways: (1) determination 
of research priorities and requirements for increasing the validity of shelter 
systems analysis, and (2) improvement of the models or formulations that are used for 
systems analysis. 

B. Procedure 

1. General - ------- 

The procedure followed in this volume is known as sensitivity analysis, 
which is a systematic determination of the effects of small changes in input 
variables. Such an analysis requires the utilization of an equation or a 
mathematical model that relates the input and output variables--such as reference 
intensity and estimates of casualties. Each input variable is varied, whilt 
all others are kept constant, so that the effects of these variations on the 
output are measured and compared. 
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Specific 


Because of the number and nature of the variables that influence estimates 
of the effects of the radioactive fallout, the establishment of a model for 
sensitivity analysis is complex. The major steps that were followed in this 
analysis are given below and are covered in more detail in Section II, which 
follows. A generalized flow chart of the procedure is shown in Figure 1, 
Research Steps in Sensitivity Analysis. 

Step I. Selection of Parameters and Values . 

Several values for each of eight factors contributing to radiological 
effects were selected as input parameters. Combinations of these values were 
used later in equations (derived while performing Step II) to provide the cases 
upon which the analysis was based. (See Section II.A below.) 

Step II. Derivation of Basic Equations to Represent Fallout Effects . 


Equations were derived for determining ERD and corresponding probability 
of fatality (P^) on the basis of the eight selected factors. (See Section 
II.B below.) 

Step III. Selection of Formulas for Determining Relative Sensitivity of 
Input Variables . 

The following relationship was established as the sensitivity index: 





Where S^ is the sensitivity index; Y is the measure of effectiveness for a 
selected set of cases; aY is an arbitrary value change in Y; Xj, is the value 
of any single input parameter in the same set of cases; and AX^ is the change 
in this input necessary to cause AY. (See Section II.C below.) 


Although it would have been desirable to relate the equations in Step II 
directly to the sensitivity ratio in Step III, it was not possible to do this 
readily without making over-simplifying assumptions. Consequently, intermediate 
steps were required. 


Step IV. Compilation of Cases for Analysis, 

Several levels of each input variable were selected and the Maximum ERD 
and probability of fatality were computed for each combination of input 


2 






























variables. Using a CDC 1604 computer and the Mainline Program which was designed 
specifically for this purpose, approximately 43,000 cases were developed. After 
those cases were excluded in which the Maximum ERD was less than lOOr or greater 
than 700r, a total of 10,143 cases were left for the sensitivity analysis. 

(See Section II.D below.) / 

Step V. Grouping the Cases and Calculation of Mean Values . 

• A tape containing the 10,143 cases was transferred to a general purpose 
data processing program written for the IBM 7072 [Reference 1] in order to group 
the selected cases into subsets for analysis. (See Section II.E below.) 

Step VI. Determination of Model Approximation Equation . 

Using the 10,143 selected cases and multiple regression, the same general 
purpose program [Reference 1] was utilized to determine a single equation 
(the principal model approximation) for the calculation of Maximum ERD as a 
function of the eight input variables. (See Section II.F below.) 

Step VII. Determination of Principal Sensitivity Indices . 

Ihe effect of each input parameter on estimating fatalities was calculated 
using the sensitivity formula (Equation 1) stated in Step III, the mean values 
of the input variables determined in Step V, and the principal model approxima¬ 
tion referred to in Step VI. (See Section II.G below.) 

Step VIII. Verification of the Principal Sensitivity Indices . 

The change in sensitivity indices was examined throughout a range of values 
using a Total Dose Model. (See Section II.H below.) These results were com¬ 
pared to those contained using the model approximation equation of Step VI. 

(See Section il.I below.) 

Step IX. Sensitivity Analysis of the Effect of Input Parameters on Shelter 
Stay Time . 

A subgroup of cases was used in an analysis of sensitivity in which the 
dependent variable was the shelter stay time required for the ERD to return 
to 80r. (See Section II.J below.) 

Step X. Evaluation of Significance of Individual Sensitivity Indices . 

The sensitivity indices of the various parameters were then evaluated in 
view of the accuracy with which these parameters could be expected to be 
established as inputs to any analysis of system effectiveness. (See Section 
II.K below.) 
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II. THE SENSITIVITY ANALYSIS 


A. Selection of Parameters and Values 

In order to obtain a representative sample of cases on which to base the analysis 
of sensitivity, a range of reasonable values was assigned to each of the input variable 
commonly used in shelter systems analysis. These values are summarized in Table I. It 
will be noted that in all cases except that for protection factor (PF), the values are 
evenly spaced. 


TABLE I 

Values and Symbols for the Input Parameters 
Used in the Sensitivity Analysis 


Parameter 

Symbol 

Values 

Reference Intensity (r/hr) 

I 

o 

300, 1500, 2700, 3900 

Time of Arrival (hrs) 

t a 

1, 4, 7 

Buildup Duration Factor 

E 

1.13, 3.64, 6.15 

Protection Factor 

PF 

2, 10, 40, 100, 500 

Time in Fallout (hrs) 

T 2 

0.0, 0.3, 0.6 

Field Decay Exponent 

Z 

-1.0, -1.2, -1.4 

ERD Recovery Fraction 

F 

0.85, 0.90, 0.95 

ERD Recovery Rate (Fraction/day) 

B 

0.020, 0.025, 0.030 


Fallout hazard (Risk) data on 200 Standard Metropolitan Statistical Areas 
(SMSA's) were used as the basis for selection of values of reference intensity and 
time of arrival. These data were extracted from an OCD source. Times of arrival 
between one hour and seven hours occurred in about 50 percent of those cases; 
reference intensities were between 300 r/hr and 3900 r/hr for 94 percent of the cases. 
(See Volume II, Figs. 31 and 32.) 

The buildup duration factor (E) is a factor by which time of ariival is multiplied 
in order to represent estimates of possible time between the beginning and cessation 
of fallout. Appendix A contains a discussion of buildup functions. 

Ihe protection factor (PF) values were based partially on National Fallout Shelter 
Survey data summarized in Reference [3]. Identified shelter spaces with PF's of 40 
and above could accommodate roughly 20 percent of the United States population. PF's 
of 40, 100, and 500 were selected to represent the three levels of PF's between 40 
and 99, between 100 and 250, and greater than 250, respectively. Factors of 10 md 2 
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were used as representatives of the PF of population in residential basements and in 
homes without basements, respectively. ; 

/ime in fallout, prior to shelter occupancy was set arbitrarily at 0.0, 0.3, 
and 0.6 hours. 

In determining equivalent residual dose (ERD), values used for the radiation 
field decay exponent and for the ERD recovery fraction and ERD recovery rate were 
varied about the generally accepted values used in systems analysis. 


B. Model Basic Equations 

The following basic equations constitute the analytical formulation of the 
model. It was upon these equations that the Mainline Computer Program (see Appendix 
B) was developed. 

1. Equivalent Residual Dose (ERD) 


ERD (n) ■ D + d 
' n r 


D n “ (1 " B) °n-l + F Cl * (t) dt 

d n * d n-l + G-F) Cl * (C) dt 


( 1 ) 

( 2 ) 

(3) 


where: 


Recoverable dose remaining during time interval n-1 to n 
(roentgens). 

III 

dose remaining during time interval n-1 to n 


of recovery for reparable damage (fraction/ 


d ■ Nonrecoverabld 
n 

(roentgens). ||| 

B ■ Fractional rat 

——-unit time). 

F * Fraction of dose which is reparable. 

(1-F) * Fraction of dose which is nonreparable. 

n ■ Number of time units after detonation, n * 1, 2, 3, ...» 
when increments (n) in time are 24-hour days. (Time is 

incremented in minutes for the first 24 hours and in days 

thereafter.) 

t * Number of hours after detonation. 

$(t) ■ Field Intensity in r/hr at time (t) after detonation. 


and: 


t(t) - 0 


for: t < T. 

** A 









where: 










sensitivity index for each input parameter. 


We begin the mathematical development by first defining the sensitivity index 
as the ratio of the percent change in the output, or dependent variable, to the per 
cent change in an input parameter [Reference 4], 

Thus, the sensitivity index (S^) of the i C ^ input parameter, measured from 
mean values, is: 

[a. gc in the dependent variable (AY)] 

c [Mean . ’ue in the dependent variable (Y)] 

5 i -- 

[Change in the l n input parameter (AK^)] 

[Mean value of the i^ input parameter (X^)] 

or: 





D. Calculation of ERD and Probability of Fatality 

The model basic equations were used with each combination of input values in 
the Mainline Program^ on the CDC 1604 computer. 

For each combination of input parameters, the time required for the buildup 
of the radiation field is calculated. Next, the dose a person would receive, if 
any, while outside the shelter and the dose at the end of the first day inside is 
computed and totaled. The dose is then split into a recoverable a^d nonrecoverable 
fraction and a test is made to determine if the probability of fatality equals 

"one." If so, the program returns to the initial state and selects another in put_ 

combi .tion. 

The ERD is accumulated and updated by recoverable and nonrecoverable fractions 
and the following tests are made on the daily increments: 

(1) If a Maximum ERD has been reached, its values and date of occurrence are 
recorded for output. 


— An additional independent variable, D,, was included to set three possible levels 
to which the ERD could return—80r, 50r, and 20r. There was a total of 42,740 
cases for all combinations of input values. 


8 - 






/ 


(2) If a Maximum ERD has been reached and the nonrecoverable ERD has not ex¬ 
ceeded the allowable dose, and the accumulated ERD is less than the 
allowable ERD, then the day is recorded for output. 

Eventually the test conditions will be met and the program will branch out of 
the day-by-day calculation of the ERD. Once out of this routine, the program will 
take the Maximum ERD value and calculate the probability of fatality and other de¬ 
pendent variables. _ _______ 

The dependent variables calculated by this program are: 

(1) Maximum ERD. 

(2) Probability of fatality. 

(3) Probability of being a noneffective. 

( 4 ) Day in which the Maximum ERD occurs. 

(5) Day in which allowable ERD occurs. 

(6) Nonrecoverable dose on that day. 

Cases were excluded in which the Maximum ERD was greater than 700r cr less than 
lOOr, thereby leaving for analysis a group of 10,143 cases (Gi.oup 0 of Table II, 
Basic and Subgroup Divisions of the Mainline Program Cases) . Thus, the sensitivity 
analysis was based on cases in which the interaction of input parameter values re¬ 
sulted in situations of maximum interest--that is, cases in which even small varia¬ 
tions in input parameters could affect both casualties and fatalities. 

E. Grouping of Cases and Calculations of Mean Values 

Subgroups of Group 0 (which contained all cases between lOOr anc 700r Maximum 
ERD) were established for each of the four input values of reference intensity. In 
addition, it was desirable to combine two of these subgroups in order to have a 
subgroup with "moderately" high reference intensity (Subgroup 20 of Table II). 

In order to examine the significance of the time-in-shelter data, another 
basic group of data was established (Group 1). Into this group were put cases in 
which the ERD would eventually reach 80r or less, because 80r is considered to be 
the threshold below which physical disability is unlikely to exist. Within this 
latter basic group, subgroups based on reference intensity, allowable ERD, and 
protection factor were established. Table II shows the basic and subgroup divisions. 

The cases were transferred from the Mainline Program by a tape to the IBM 7072 
computer. A general purpose data processing program (TSAR) at Duke University was 
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TABLE II 


Baste and Subgroup Divisions of the Mainline Program Cases 


BASIC CROUPS SUB-GROUPS 

NUMBER 

SORTING PARAMETER 



OF CASES 


Group 0 


10,143 

Reference Intensity (r/hr) 

(parameter Combinations 
with : 

5 

lOOr < Max. ERD < 700r) 6 


2,115 

3,108 

I - 300 

1° - 1500 

7 


3,431 

1° - 2700 

8 

’ ..."... 20 

of 

(sum 

7 & 8) 

1,48? 

4,920 

1° - 3900 

1° - 3063 
o 

Group 1 


3,481 


(Cases which return to 
an allowable ERD < 80r) 



Reference Intensity (r/hr) 

9 


909 

I - 300 

10 


1,041 

1° » 1500 

11 


1,058 

1° - 2700 

12 


473 

1° - 3900 

0 




Allowable ERD (r) 

2 


2,069 

D. - 80 

3 


1,204 

d£ - 50 

4 


208 

d£ - 20 




Protection 

13 


987 

PF * 2 

14 


1,291 

PF - 10 

15 


795 

PF * 40 

16 


321 

PF - 100 

17 


87 

PF - 500 




Multiple Sorting Parameters 

18 


352 

D 4 - 80r 




0 < P 2 < .5 
t 2 - 0 

19 


240 

D. - 50 r 

4 


0 C p 2 < .i 

t 2 - 0 
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used for the grouping of data into sets 'RofJrence l . 

The same program was used for determining the mean values of the input parameters 


for the basic groups of cases and for each of] 
of the Basic Parameters , shows the pertinent 
groups. 


the subgroups. Table III, Mean Values 
heans for each of these groups and sub- 


F. Determination of the Principal Model Approximation 


For purposes of computation it was desir 


able to have a single equation that 
approximates the Mainline Program computer model. Using a stepwise regression to 
obtain this, three basic functions (quadratic 


, logarithmic, and exponential) were 
fitted to the sets of data produced by that mlodel. This was accomplished also by 


the TSAR general purpose program. Each funct 


ion was fitted by th 2 step procedure 


*[Reference 11. Iflquation 5 best fitted the data and was made the principal model 
2 / ' 

approximation :*- 


f (X) => Max. ERD ■ 
(-2.7468) 


(.932 + .0001499PF 


- .0446F + 5.354B - .006493E) 


T (.4346 - .7954Z) pp (-1.0949 + .04181T 2 4- .06268Z) 
z (-3.074 - .1161E) b (- 1.1553) j 


(5) 


Th 


e statistical characteristics of the principal model approximation equation 


(1) Multiple correlation coefficient: 0.9991. 

(2) Standard error of the estimate: 0.02421 (i.e., c of ERD + 6r). 

. ' ! .. 

(3) Value of F-ratio: 886.785 (a further indication of the high significance 

I 

of the equation). 

(4) Data Group 0 (Table II). 

G. Calculation of the Sensitivity Indices 

We now have a formula (Equation 5) relating the input variables to the dependent 
variable Maximum ERD, and a formula for the sensitivity index (Equation 4). To 
calculate the sensitivity indices at the mearjs, Equation 4 is placed in differential 
form: 


—I Identification of the variables is in Table I. 
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Individuals in 500 PF shelters only reached or exceeded the Maximum ERD if they receive appreciable dose outside. 
Their time-outside average was twice any other group. 






£1 

Y 


(^) 


£X X (Y)" 1 

&x, x i w 


Lee: 


Y - f(X) 


s i -‘aH < £ p » _1 


( 6 ) 


Note that f(X) is the principal model approximation (See Equatioti 5): 


^1 ^2 ^3 ^5 

f(X) • K I T. PF Z B 


(7) 


where: 


C 1 .932 + .0001499PF - .0446F + 5.354B - .006493E 
C 2 - .4346 - .7954Z 
C 3 - -1.0949 + .04181T 2 + .06268Z 
C 4 - -3.074 - .1161E 

C 5 * -1.1553 
-2.7468 


Ihe partial derivatives of Equation 5 with respect to each independent variable 
(X t ) are: 

mc/vN C.-l C, C- C, C, 

»KC. I 1 T. 2 PF 3 Z 4 B 5 
c*t 1 o A 

o 


- 3 ' 1 ‘ C4 » s 

- K K I Cj T ° 2 PF C3 Z C4 B C5 
5T 2 K K 3 o T A 

C, C„ C C,.-l C c 
PF J Z 


aaa - K K, r 1 t / 2 pf “ 3 z ~ 4 * B ~ 5 


>Z 


1 4 o 


c, C„ C, C, C- 

21^- * K K. I 1 T 2 PF 3 Z 4 B 5 
BF 5 o A 
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F - .9 
B - .025 
E - 2.57 

and using the Equation 4, we calculate the Sensitivity Indices (S^). 


TABLE IV 

Principal Sensitivity Indices 


Index 

Variable Denoted by Index Subscript 

Value 

s z 

Field Decay Exponent 

-4.22 

S PF 

Protection Factor 

-0.96 

S I 

o 

Reference Intensity 

+1.02 

St a 

Time of Arrival 

-0.50 

S F 

ERD Recovery Fraction 

-0.30 

S B 

ERD Recovery Rate 

-0.14 

S E 

Buildup Duration Factor 

-0.18 

S T 

l 2 

Time in Fallout 

40.05 


It is upon these values that the conclusions and recommendations will be 
developed. 


The above sensitivity indices are obtained for Maximum ERD as the dependent 
variable. However, over the dose range where probability of casualty is neither 
0 nor 1, che relation between probability sensitivity (PSp and dose sensitivity 
(S^) Is given by the relation: 

PS i " S i (1 * p> <»> 

where p ■ probability of casualty = a • Dose + b. Hence, the relative rank of the 
S f remains the same when evaluated at a given set of input variables. 

H. Total Dose Model 

The major computational difficulties in determining the above sensitivity 
indices arise because of the use of equivalent residual dose instead of total dose 
as a measure of biological effect. This section will describe a simpler model which 
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allows the calculation of the sensitivity indices analytically. In this model the 
dose rate buildup function is assumed to be linear, and total dose is used as the 
dependent parameter. 


A new parameter, T o , is introduced in this approach. T q is the time period from 


entrance into shelter until the cutoff of the total dose calculations. Ihus, 

(T^ + T 2 ■]+ T q ) represents the complete time period from time of detonation until com¬ 


pletion of the total dose calculation. Using the parameters previously defined (see 


Table I) and the new parameter T^, the total dose D is given by the expression: 


D ‘P? 


2T Z+1 <1+E) Z E 


2 2 2 
(PF-1) T, + E T. 


<1 




(T a + T 2 + T o ) UZ - T a 1_Z (1+E) 


"1 


The sensitivity indices can then be calculated directly 
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The total dose sensitivity indices (S^) can be readily converted to probability 
of casualty sensitivity indices (PS^, if desired. In the region of interest 
(probability neither zero nor one), the dose is related to probability of casualty 
(p) by a linear relation: 

p * aD + b. (18) 


ps i * lx- i ~ h « • ?> 


Only dose sensitivity will be discussed in the remainder of this analysis, since the 
relative rank o' both S and PS are the same for a given set of values of the inde¬ 
pendent variables. 

In Table V, Comparison of Sensitivity Indices , the values of the S^ obtained 
from Equations 7 through 13 by evaluation at the mean values of the variables listed 
in section G above are compared to the values of S^ listed in Table IV (calculated 
from the ERD model). 
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TABLE V 

Comparison of Sensitivity Indices 


Index Variable Denoted by Subscript ERD Model Value 


S_ Field Decay Exponent -4.22 


otal Dose Model Value 


Field Decay Exponent 

-4, 

.22 

-4. 

,36 

Protection Factor 

-0, 

.96 

-0. 

,99 

Reference Intensity 

+1, 

.02 

+1, 

.00 

Time of Arrival 

-0, 

.50 

-0. 

50 

! 

Buildup Duration Factor_ 

_-0. 

.18 

-0. 

,25 

Time in Fallout i , 

+0. 

.05 

40. 

.02 


In Figures 2 thru 6, the values of the sensitivity indices are plotted as a 
function of the indexed variable. While the. rest of the independent variables are 
held constant at the mean values given previously, two values for PF (39.1 and 500) 
and two values of T q (96 and 240 hours) are plotted for each curve. Figure 7 shows 
the increase in the sensitivity index for (time in fallout) as Z (the decay rate) 
increases. Figure 8 shows the dose sensitivity to T as a function of T for the 

OiO 

two protection factors. This curve indicates that total dose jis not particularly 
sensitive to T q when T q is measured after 100-200 hours, especially for high PF 
shelters. In fact, where Maximum ERD occurs within the first jfew days, it can be 
expected that the results from the total dose sensitivity analysis model will be 
essentially equivalent to the results from the ERD model. 

I. Validity of Ranking the Sensitivity Indices 

In Tables IV and V, the sensitivity indices were evaluated for one set of input 

i 

variables. However, the relative values of these Sj, might change when evaluated for 
radiological environments other than the one defined by the mean values in section G. 
To investigate this problem, the following conclusions were drawn based on Figures 
2 through 8 and other calculations of sensitivity indices: ! 


(1) Sj , sensitivity index for reference intensity (I 0 )> jis 1 over the entire 

range of values; thus, I q remains a relatively important parameter in all 
cases. j 

(2) S n _ is essentially constant at 1 for most cases of interest; thus, protec- 

PF 

tion factor (PF) is a relatively important parameter in all cases. (See 
Figure 2). 
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Fig. 2. Sensitivity Index of PF as a Function of PF 



























































Fig. 4. Sensitivity Index (S_) as a Function of Buildup Factor (E) 








































Sensitivity Index (S„, ) as a Function of Time in Fallout 
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Fig. 6. Sensitivity Index (S_) as a Function of Decay Constant (Z) 
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(3) S,, remains high for all cases, thus Z is extremely critical in vulnerability 
analyses. (See Figure 6). 

(A) For analysis observing only Maximum ERD reached while in primary shelter, 

ERD concepts and a total dose approximation of Maximum ERD yield essentially 
the same results in sensitivity analysis. (This does not imply that they 
yield the same results in casualty calculations.) 

(5) can become quite large for short times of arrival of fallout (T^) and 
high protection factors. In areas where this combination of environments 
exists, T^ must be treated carefully in vulnerability analyses. (See 
Figure 3). 

(6) In most cases of interest, the other sensitivity indices remain appreciably 
less than 1, though some change significantly as the parameter values are 
varied. 

Analysis of Time Required in Shelter 
1. Procedure 

In the preceding analysis, the independent variable was Maximum ERD. It 
is also important to know how the input variables may affect the time required 
in shelter after the beginning of fallout. 

In the analysis of time required in shelter, only those cases from the 
Mainline Program were considered in which the probability of a noneffective 
(P 2 ) was greater than zero and less than 0.5. Thus, the largest Maximum ERD 
to be considered was 200r and the smallest Maximum ERD was lOOr. (See the dose 
response relationships in Section II.B.2.) Another condition was that the ERD 
returned to 80r. It was assumed that the reduction of ERD to this level would 
mark approximately the end of the continuous stay in shelters. Time in fallout 
outside shelter was also excluded from consideration. 

The TSAR general purpose program was used to set up subgroups 18 and 19 
(see Table II) and to develop an equation by multiple regression. The equation 
given below determines the number of days in shelter until ERD returns to 80r. 

When: lOOr £ Maximum ERD < 200r, 








( 20 ) 


* t I 

f(X) - N 2 «* 4.47 - 44Z + 2.26 p + 1.6D, - I.55(Z x ~) 

+ 926(Z x B) * .7920(D^ x F> - 19.8(D 3 x B) 

if 

Multiple Correlation Coefficient: 0.9717 
Value of F-ratio: 830 

Where: 

N 2 ■ Duration of shelter phase until ERD falls to 80r 
Z ■ Field decay exponent 

J o ■ Reference intensity/protection factor 
PF 

«* Maximum ERD prior to N 2 

B ■ ERD recovery rate 

F ■ ERD recovery fraction 

The mean value of stay time for the cases considered (Y in the sensitivity 
formula) was found to be 47 days. 

The sensitivity indices were calculated using Equation 20, the differentia, 
form of Equation 4; 


& 



c , MOO i 

i ?X i ’ 47 

witn: 

N 2 * Y * f(X) (see Equation 20), 

N 2 * Y = 47 days. 


( 21 ) 


2. Results 

The sensitivity indices for the input parameters are shown in Table VI, 
Sensitivity Analysis of Time Required in Shelter . It is seen that changes in 
the decay exponent (Z) are the most significant in affecting the length of stay 
in shelter. 
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TABLE VI 


// 

/ 
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Sensitivity Analysis of Time Required In Shelter 


Parameter 

1 

Mean Value 

*i 

Gradient 

£f(X) 

dx t 

Sensitivity Index 
(Equation 21) 

S i 

Z 

1.260 

- 144.000 

-3.860 

F 

0.914 

- 144.000 

-2.220 

°3 

144.000 

0.373 

1.140 

B 

0.025 

-1700.000 

-0.918 

RP 

79.600 

0.331 

0.560 


The independent variables (which affect the length of stay in shelter) are 
ranked on the sensitivity index in decreasing importance as follows: 


1. Field decay exponent, Z “3.9 

2. ERD recovery fraction, F » 2.2 

3. Maximum ERD, * 1.1 

4. ERD recovery rate, B » 0.9 


5. Reference intensity/protection factor ratio, RP *• 0.6 

Of equal interest to the parameter sensitivity analysis is the actual 
timing of Maximum ERD and the number of days' stay in the primary shelter before 
recovery to an allowable dose is accomplished. For a shelter PF of 73, the 

Maximum ERD of 144r occurred on the 7th day after the attack. On the average, 

3/ 

47 days elapsed before the ERD recovered to the allowable dose level of 80r.— 
Figure 9 shows the frequency of occurrence plotted against the duration in days 
before reaching an ERD of 80r. It can be seen that there are cases in which 
stay time would exceed 100 days if the ERD of shelter occupants is to return 
to 80r. 

K. Evaluation of Significance of Sensitivity Indices in Shelter Systems Analysis 
1. General 

The objective of shelter systems analysis is to predict as accurately as 
possible the adequacy of existing or potential protective measures against 

3/ 

— If instead, an allowable level of 50r (Subgroup 19) were chosen, then the average 
length of stay would increase by 33 days (from 47 to 80). 
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Days for Return of ERD to 80r (N 2 ) 

FLg. 9. Frequency of Occurrence of Stay Times vs Days for Return of ERD 
(for Cases in Which Max ERD < 200) 










fallout. In general, the better the parameters that determine the adequacy are 
known, the better the predictions will be. 

The sensitivity indices derived in the above analysis give an indication of 
the effects of not knowing the various input parameters exactly; the analysis 
does not yet take into account how accurately these parameters might be known. 

In the following sections, an attempt will be made to identify the approxi¬ 
mate magnitude of uncertainty to be expected in the eight parameters considered. 

2. Reference Intensity and Time of Ar-ival 

Of the eight parameters that were considered in the sensitivity analysis, 
reference intensity and time of arrival are largely dependent upon the choice 
of strategy in the attack; consequently, these parameters cannot be known 
accurately. Desired ground zeros and the occurrence of air or surface burst are 
the choice of the attacker. 

At any given point, the reference intensity and time of arrival are also 
affected by wind speed and direction. 

The RISK type programs are attempts to represent the combined effects of 
strategic decisions and wind conditions. In the choice of input variables for 
the sensitivity analysis (section II.A), it was seen that reference intensities 
between 300 and 3900r/hr and times of arrival between 1 and 7 hours seem to 
represent the range of intensities with which fallout-only shelter systems are 
primarily concerned. Variations from the means of these input parameters of 
+ 75% would seem to be a representative estimate of the probeble range of values 
for reference intensity and time of arrival in an analysis of a shelter system. 

3. Buildup Factor 

The values of the buildup multiplying factor used in the sensitivity 
analysis ranged from 1.13 to 6.15. These extreme values vary from their mean 
by about 707. and this is proposed as a range of values that could be expected 
in systems analysis. 

4. Protection Factor s 

The National Fallout Shelter Survey (NFSS) produced data on the present 
number and the protection factors of shelter spaces. However, there are various 
reasons why these data alone may not represent the actual shelter posture of the 
population. 
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First, based on a limited sample of shelter [Reference 5], it appears that 
procedures used in the Phase 1 shelter survey generally underestimated the 
actual protection factor by a considerable amount. The current updating of 
NFSS data is gradually improving the accuracy of information on shelter status, 
but a complete reevaluation of all NFSS shelters is not contemplated by OCD. 

At the present time the average potential protection factor of the United States 
population may be considerably more than presently represented in Phase 1 or 
Phase 2 data. However, ever, if the protect tur. factors w\.re accurately deter¬ 
mined, there still would remain uncertainties as to shelter utilization that 
would influence the results of shelter analysis. The question of shelter 
utilization was presented in Volume II where it was considered in terms of 
actual SMSA's. (See Volume II, A Sensitivity Analysis of Selected Parameters 
Based on 8 SMSA's.l 

Consequently, random variations of 807. about the mean can be considered 
reasonable in estimating the significance of the sensitivity index for protec¬ 
tion factor. 

5. Time in Fallout 

Estimates of length of time outside in fallout are arbitrary. By assuming 
adequate warning, no one would need to be outside. However, there may also be 
cases where large numbers of the population are in fallout for periods ex¬ 
ceeding the 0.6 hours arbitrarily used in the sensitivity analysis. 


Although the sensitivity index for time in fallout is small, the importance 
of this parameter in systems analysis must be determined in conjunction with 
actual locations and categories of shelters. The combined effects of these are 
analyzed in Volume II. 

6. Decay Rate - .-.—----—. .- 

Ihe decay rate may vary with the construction of the weapons and to a lesser 
degree with the nature of the surface materials at ground zero. Variations 
probably will not be more than plus or minus 25 percent of the 1.2 rate of 
decay, which is based on empirical data. 

7. Biological Factors 

The fraction of the dose that is reparable (F) and the rate of recovery 
(B) are important parts of the Blair Formula for calculating ERD. However 
according to Reference [6], a firm experimental basis appears to be lacking 
for the values presently used. The author in the same reference also questions 








whether the constants are independent of time. 

In this study, the average value of F, the recovery fraction, is taken as 
0.90; variations of 10 percent would seem likely. Ihe average value of B, the 
recovery rate, is 0.025 and probably could vary by 70 percent. 

8. Comparison 

We then may compare the sensitivity indices with the rough estimates of the 
accuracy to be expected in input variables. These are sunmarized in Table VII. 
An approximation of the relative importance of each parameter in contributing to 
total uncertainty in shelter analysis is obtained by multiplying the range of 
uncertainty by the sensitivity indices. 

TABLE VII 

t- 

Comparison of Sensitivity Indices and Uncertainty in Input Parameters 


Approximate 

S (ERD S (Total Relative 

Variable Model) Dose Model) Uncertainty Importance 


Field Decay Exponent 
Refererce Intensity 
Protection Factor 

Time of Arrival 
Buildup Duration Factor 
ERD Recovery Rate 

Time Outside in Fallout 
ERD Recovery Fraction 


-4.22 

-4.36 

+1.02 

+1.00 

-0.96 

-0.99 

-0.50 

-0.50 

-0.18 

-0.25 

-0.14 

— 

+0.05 

+0.02 

-0.30 

— 


25% 

1.09 

75% 

0.75 

80% 

0.80 

75% 

0.35 

70% 

0.18 

70% 

0.10 

100% 

0.05 

10% 

0.03 


lift. CONCLUSIONS AND RECOMMENDATIONS 


A. Conclusions 


Considering both the sensitivity of each parameter and the range of uncertainty 
that each is likely to have, it is possible to draw some conclusions as to their 
effect on the results of systems analysis. 

In analyses for estimates of Maximum ERD, the following have the most signi¬ 
ficance: Field decay exponent (Z), reference intensity (I 0 )> and protection factor 
(PF). Having less significance are time of arrival (T^) , ERD recovery fraction (F), 
ERD recovery rate (B), and buildup duration factor (E). 
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Time in fallout (T£) has a low sensitivity index. However, extending the range 
of values that were assumed for it beyond 0.6 hours may have a considerable effect 
on the magnitude of this number. The importance of time in fallout was better shown 
in Volume II, where it was analyzed as a factor in shelter utilization. 

In analysis of required time in shelter, field decay exponent (2) and ERD re¬ 
covery rate (B) are the most significant. 

Of the above parameters, reference intensity and time of arrival are almost 
entirely dependent upon the strategic choice of the attacker and wind conditions. 
Consequently, better estimates of their effects in systems analysis is dependent 
upon military intelligence more than upon additional research and analysis. 

The importance of the field decay exponent in determining both fatalities and 
shelter stay makes it important to have continuing analysis of the validity of the 
1.2 value. There may be appreciable gains in accuracy possible since the exponent 
is primarily dependent upon the physics of the radioactive decay process. 

Increased knowledge of the shelter protection factors available to the population 
is important for the making of valid systems analyses. The effect of underestimates 
of protection factors can be taken into account by using a multiplying factor. 

However, it is important to define the factor more closely than is possible now. 

This would not necessarily require large expendir.res for obtaining new data or for 
correcting data that now exist. Statistical analysis based on the parts of the 
survey that have been updated or examined in detail could provide a much better 
estimate than evidently now exists of the overall shelter that would be available to 
the population. 

Identification of the biological recovery process is particularly important in 
estimating time required in shelter. Possible approaches to reducing tne range of 
uncertainty of this input are given in Reference [6], 

Finally, efforts to define more closely the limits of any input variables should 
be viewed in the light of their relative effect on output data. For example, attempts 
to define protection factors accurately may make a change in the relative importance 
of that parameter that is small compared to the relative importance of reference 
intensity, which cannot be accurately predicted by any research means. 

The sensitivity index of the buildup duration factor (E) is based on "moderate" 
conditions where the time of arrival approximates 3.3 hours and (E) approximates 2.6. 

In observing the value of the sensitivity index S^,, it was found that the index is 
quite dependent on changes in (E). Some further analysis of the effect of the build¬ 
up duration factor is discussed in Appendix A. 
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The analysis based upon return of ERD to 80r (the threshold of appearance of 
clinically detectable symptoms of radiation sickness), showed that shelter stay times 
of much more than two weeks could be necessary. It should be noted that this analysis 
included only cases in which the Maximum ERD was 200r or less. Without these restric¬ 
tions, even longer required stay times would result. This is very significant in 
planning for the stocking of shelters, recovery planning, and for the establishment 
of policies for postattack management of supplies. 

In summary the conclusions from the sensitivity analysis are: 

(1) Casualty calculations are quite sensitive to errors in the field decay 
exponent. The sensitivity of dose or casualties to the field decay rate 
remains at approximately -4.0 over the range examined. (The decay rate 
was varied from Z «* 1.1 to Z •* 1.4.) 

(2) Sensitivity to variations in fallout reference intensity and protection 
factor are high.over the whole range of parameter values. Hence, precise 
knowledge of the fallout shelter posture and the fallout reference inten¬ 
sity is much more essential to accurate vulnerability analysis than the 
remaining parameters in most cases of interest. 

(3) Sensitivity to time of arrival of fallout can be quite high, in some 
radiological environments. 

(4) Sensitivity of casualty computations to the remaining parameters is low 
in most cases of interest. 

(5) Expansion of the sensitivity analysis to include parameters other than 
fallout, which define the total casualties from a given attack on the 
United States, is necessary before further conclusions concerning a 
national vulnerability analysis can be drawn. 

B. Recommendations 

1. Vulnerability analyses should employ protection factors computed by best 
available methods, and research and/or surveys to improve protection factor data 
should be encouraged. 

2. Because oi the sensitivity of systems analysis results to the field decay 

-1 2 

exponent, continuing analysis of the validity of the t ’ decay law should be made. 

3. Additional study is required to establish the sensitivity of fatalities, 
casualties, and dose to duration of shelter stay. 
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4. Sensitivity analysis should be extended to include the parameters defining 
the effects of blast and alternative measures of effectiveness (casualties by type, 
lose distribution of survivors, etc.)* 

5. Sensitivity analysis should be applied to identification of the important 
cost/effectiveness parameters used in the budget allocation model (see Volume I, 

k Cost/Effectiveness Computer Procedure for Optimum Allocation of Fallout Shelter 
System Funds Under Uniform or Variable Risk Assumptions) . 
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Appendix A 


The Intensity Buildup of a Radioactive Fallout Field 


In calculating the radiatior dose received at a point in the fallout field, it 
is necessary to determine the rate at which the fallout particles arc deposited from 
the 1 cloud. The easiest method o i representing this buildup is to plot the fraction 
of jtotal fallout mass deposited £,s a function of time. If is the time the first 


particles arrive, and is the tjime the last particles are deposited, then the 
duration of the buildup period, - T^,is approximated by a relation of the form 


where E is the buildup duration factor referred to in this volume. 

I: _ 

Figure A-l presents the buildup function ''Reference A-l" which was obtained 
from weapons tests. In addition^ two approximations to this function are plotted 
in ’Figure A-l for = 1 and E =2.5. The analytic form for these approximate 
curves is given by Equations A-2 and A-3. 


‘ V 

F <t)j- ? - cos ~) 


Equation A-3 gives a somewhat better fit of the curve from Reference r A-l', 

hence this buildup function was used in the computer model that calculated Maximum 

i 

ERD. However, the tctal dose expressions cannot be integrated in closed form if 
Equation (A-3) is used, hence Equation (A-2) was used as the buildup function for the 
analytical total dose model. AsJ might be expected from Figure A-l, little difference 
was found between sensitivity indices from, the two models. 
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[Afl] Miller, C. F. Fallout a nd Rad iological Countermeasur es, SRI Project IM-4021, 
Vol. I, January 1963. 
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Fraction of Total Fallout Mass Deposited 



Time after Detonation 
Fig. A-l. Fallout Mass Buildup Functions 
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Appendix B 


The Mainline Model and Computer Program 


I. THE OPERATIONAL TRANSFORM OR MAINLINE MODEL 


A. Attack Environment Formulation 


1. Introduction 

Ihe first step in formulating the Mainline Program is proper characteriza¬ 
tion of the attack environment. Ihe interrelation between the initial intensity, 
field buildup, the field decay, and the field attenuation by shelter, are the 
four functions that predict the dose rate. To determine their relative im¬ 
portance, it becomes mandatory that tnis interrelation is accurately formulated. 

The underlying assumptions throughout this section are: 

(a) Ihe Reference Intensity of the radioactive field (in roentgens per 
hour) will conform to the measurements given in the NREC RISK-type 
data (i.e., an intensity (I q ) at H + 1 hours) as well as the distribu¬ 
tion of intensity and arrival of fallout obtained also from NREC RISK- 
type data. 

(b) Ihe measure of protection obtained by the emploj'ment of shelters will 

be the Protection Factor (PF) obtained from "The Fallout Shelter Surveys" 
_Phases 1 and 2._■_____! _ ....... 

2. Decay of Field Intensity 

Since a radioactive particle decays with time, our problem is to determine 
the appropriate mathematical representation, when applied to residual fields 
used in Civil Defense Systems Analysis. Ihe widely used function shown in 
Figure B-l and Equation (B-l) will be used. 


9(t) - I 0 t 
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where 


0(t) ■ Field Decay Function, in r/hr. 
I Q « Reference Intensity, in r/hr. 
t ■ Time in hours. 

Z • Empirical Constant. 



Fig. B-l. Decay of Radioactive Field 

The two important constants introduced with this function are: the 
empirically derived exponent (Z), and the reference intensity (I q ). The refer¬ 
ence intensity is one of the two interface parameters connecting the fallout 
shelter analysis model with the '*war games" which establish the probable attack 
environment. 


3. Buildup of Field Intensity 
a. Introduction 

The problem now is to characterize analytically the field intensity 
during the buildup phase. The fact that there is no existing data within 
our files on which the usual methods of curve fitting can be applied, forced 
us to use a more general form of analysis. Our analysis uses what little 
significant information there is and employs as much intuition as seems 
necessary, useful, or feasible. Although there are two relevant functions 



(B-2) 


and 


t 

e 


5T. 


7 


(B-3) 
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they are incompatible over a wide range of times of arrival. The relation 
time of peak intensity equals twice the time of arrival--is taken as being 
more nearly valid. In the field buildup, two distinct phenomena are 
present. The first is the arrival (at ground level) of the radioactive 
particles. The second is the independent phenomena of the radioactive 
decay of these particles. Thus, characterizations of these two phenomena 
will be independent, and the field intensity represented by the product. 

The hypothesis implies a particle buildup function of the following 
general form given in Figure B-2, which, when integrated, yields the trine 
dependent arrival function given in Figure B-3. 



Time t -► 

Fig. B-2. Particle Arrival Rate 



Fig. B-3. Accumulation of Particles 
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The produce of the function* represented in Figure B-l and Figure B-3 
would give the field intensity shown in Figure B-4. 



Fig. B-4. Field Intensity vs. Time 

b. Formulation 

From several general forms a displaced Cosine function was chosen 

for simplicity and compatibility with the function (t ■ 2T.) to represent 
& P A 

the fallout particles accumulated (Reference Equation B-4). 


F(t) 


*[’- 


cos (- 


*) 


(B-4) 


where 


and 


T, < t < T 

A - - C 


1 ■ T C • V 

-z 

The field decay function (i.e., t ) times the buildup function times the 

reference intensity yields the final representation of the field intensity 

over time. For time between T. and T : 

A C 


KO - I 


si |\. „>J 


(B-5) 
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Thus, the characterization of the fi< Id intensity is in two parts; the 
first during accunulation--fron T (time of arrival) to T c (time of 
cessation). The second, after all fallout has occurred (Reference Equation 
B-l), describes the field. 

In order to comply with T = 2T,, the value of must be fixed. The 
p A 

following are values of for three values of Z: 

Z - 1.0 t) - 1.133T a 

Z - 1.2 tv - 1.175T. 


4. Attenuation of Field Intensity 

a. The Problem 

The population or an individual may move through a series of shelter 
conditions over time. The duration of each shelter phase, although 
nominally under operational control, has well-defined bounds based upon 
the attack environment. The initiation of the recovery phase depends upon 
minimizing the stay in shelter. Thus, the influence of shelters when 
modeled must take into account the operational requirement to minimize 
shelter stay time. 

b. Formulation 

Fallout shelters are classified by a protection factor (PF). The 
protection factor function operates on the unprotected radiation field 
to yield the attenuated field that exists within that shelter. 

There are four major operational phases (Reference Figure B-5). 

These are: 

(1) An initial unsheltered condition which corresponds with a person 
moving through a fallout field before entering shelter. This 
phase will be called "time outside.It is related to the 
efficiency of warning systems, drills and training, and shelter 
assignment plans. 


(2) Phase 2 is the primary shelter phase. In most current models it 
is the only one considered. Operational plans usually specify 


—■ The symbol used in this report for initial time outside in a fallout field is 
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Field Intensity (r/hr) 


a two-week stay, and they are itocked on that basis. The control¬ 
lable variable is the length of stay in the phase. Our model 
considers the length of stay dependent upon an objective of the 
third phase and is subject to explicit calculation. 

(3) Phase 3 is operationally associated with emerging from shelter to 


begin recovery. The objective 


(4) Tne fourth phase is the return 


in this period is to obtain produc¬ 


tive effort (usually at a PF oJ r 1) without any physiological 
radiation symptoms. To do thii introduces a decision parameter 
(an allowable ERD), an operational parameter, and the equivalent 
PF required. The allowable ERD lies between the limits of 
clinical effects (25r) and symptomatic effects (lOOr). 


to normal. It begins when the 


individual's recovery rate equals the dose received in the un¬ 
sheltered field. It is characterized by the lack of control 
variables, decisions required, etc. In other words, the effects 
of normal biological recovery processes exceed the effects of the 
radiation fields. j 

The problem domain investigated in this volume extends only to Phase 3, 
time (T^). (Practical considerations precluded extensive investigation of 
Phases 3 and 4 under the present contract.) 



Fig. B-5. Attenuation of Field Intensity 
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In Fig. B-5, the solid line represents the sheltered point intensity. 
In Fig. B-5, the dotted line represents the field intensity if: 

PF^ « PF 2 ■ PFj " 1. 

where 

Phase 1 PFj(t) * 1, 0 < t < T 2 

Phase 2 PF 2 (t) - PF, T 2 < t < T 3 

Phase 3 PF 3 (t) ■ Equivalent PF(t), T 3 < t < 

Phase 4 PF 4 (t) - 1, < t. ' . 


Attack Environment Model 


The attack environment radioactive field intensity at any time (t) is the 
product of: the Reference Intensity (I Q ); the particle buildup function (F(t)); 
the field decay function (-"(t)) ; and the protection factor function (PF(t)). 


«(t) - R.. P(t) • (B-6) 

Or by substitution: 

;<t > ■ ?f; V [* - cos ( ~r ”’] “ he " < e £ t a + -' t a 

*(t) - 0 when t < (B-7) 

l Q - 

5(t) - t when t > T A + rjT^ 

where: 

(t) * Field Intensity in r/hr at time t after the originating 
explosion event. 

I o * The reference intensity in r/hr, one hour after the explosion 
event. 

t ■ Time in hours measured from H. 

T. - Time of arrival measured from H hour. 

A 

PF, * Protection factor at time t. 
k 

Z - empirically derived field decay constant. 

n ■ T„ - T. ■ fallout interval. 

1 C A 
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Biological Response 


1. Introduction 

The state, efficiency, and outcome of a shelter system analysis is measured 
by the biological impact, the usual criteria being the number of fatalities or 
lives saved. Fatalities are calculated by using a probability of response to 
the total radiation dose received. 

At the present state of understanding, the equivalent residual dose model 
is the best available formulation of biological response to dose received. 

2. Equivalent Residual Pose 

To represent analytically the absorbed gamma radiation field in terms of 
a time dependent dose "state" requires taking into account the observed phenomena 
of biological recovery. The ERD model relates the field intensity and biological 
recovery rate to give a single "dose" value, approximating the physiological 
state of a human. 

The ERD at any point in time consists of two parts. The nonrecoverable 
portion is the accumulative amount of the nonrecoverable dose received over 
time. The recoverable portion has two additive parts; one is a fraction of the 
field contribution within the unit time considered, and the second is the 
recoverable dose carried over from the previous period. The basic statement 

of ERD (recoverable dose (D ) in time increment, n) is the sum of that not yet 

“ th 

received from previous periods and the new contribution during the n time 

period. 

D n - (l-B)D n _ 1 + F / n «(t)dt (B-8) 

n-1 

where •__;____ 

B is the percent per unit dose recovered in period n, and F is the 
recoverable fraction. 

Nonrecoverable dose (d ) in time increment (n) is the accumulative sum of all 
n 

nonreparable doses: 

d n " d n-l + (1 ' F) /" ®(0«lt. (B-9) 

n-1 

The practical problem of describing the function of the ERD(n) in a closed, 
well-behaved form requires that it be approximated by finite differences. 

Finite difference problems, which require recursive solutions, lead to computer 
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Implementation. A solution is to use the basic difference equations with n 

in one day intervals and with d being computed daily, then added for the 

^ th 

equivalent residual dose on the n fc day. 
where 

D ■ Recoverable dose in time interval n-1 to n, 

d ■ Nonrecoverable dose in time interval n-1 to n, 
n * 

I « Reference intensity, 

F “ Fraction of dose which does reparable damage, 

(1-F) ■ Fraction of dose which does nonreparable damage, 

B - Daily fraction of dose repaired, 

T. ■ Hie time of arrival of the first measurable radioactive fallout, 

A 

Z * Decay constant, 

PF * Protection factor. 

3* Probability 

The problem is now to relate the equivalent residual dose which an individ 
ual accumulates to the probable biological state. The usual specifications for 
radiological induced states are fatalities and casulaties, with fatality being 
a type of casualty. Because this causes some confusion, we will introduce the 
term noneffective. For a given Maximum ERD, we will then have three probabili¬ 
ties: 

(1) Death. 

(2) Noneffective--The condition of being alive but physically unable to 
perform a productive task. 

(3) Normal and Marginal Effective—Individuals with acceptably low Overt 
physiological symptoms (below lOOr). 

The following linear functions are used to approximate the empirically derived 
probability of casualty vs. ERD curves (reference Figure B-6). 

a. Probability of Fatality 

Pj - 0 when ERD < 200r 

Pj _ ■ QfeXj lll ~| gg ll 2 00 when 200r < ERD < 700r 

Pj - 1. when 700r < ERD. 
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2 100r 200r 300r 400r 500r 600r 700r 

0 * 


Equivalent Residual Dose (ERD) 
Fig. B-6. Dose Response 



Probability of Noneffectiveness 



when ERL < lOOr 
when lOOr < ERD < 200r 
when 200r < ERD < 300r 
when 300r < ERD < 700r 
when 700r < ERD. 



(B-U) 


(B-12) 


Mainline Program 


1. introduction 


The models, previously described, are combined into a single nK>del--a 
computer program. The finite difference method, particularly adaptable to 
computer progTamning, is necessary because of several functional discontinuities 
which preclude integrating the individual models simultaneously. In translating 
algebraic formulas into FORTRAN coding, a change in notation is required. To 
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assist in identifying corresponding U ra 1 , a Glossary of Syr.brls use-1 in the 
programming, and the corresponding symbols found in mathematical formulation 
are listed in Table B-I. 


TABLE B-I 
Program Symbols 


Parameters 

Model Symbols 

■Fro granting Symbols 

Reference Intensity 

V 

R 

Time of Arrival 

. T A 

I/. 

Time Outside 

T 2 

. 

Allowable Dose 

D 4 

!>- 

Field Decay Constant 

Z 

Z 

Recovery Fraction 

F 


Recovery Rate 

ft 

0 

Buildup Duration 


f:;a 

Protection Factor 

PF 

PF 

Protection Factor Categories 


K 

Maximum No. PF Categories 


KMAX 

Time (the variable of 
integration) 

Hours 

Days 

n 

T 

Total (ERD) Dose 

ERD <n) 

D 

Recoverable D.-se 

D 

n 

D1 

Konrecover<iMe Dose 

ci 

n 

D2 

Day Maximm ERD Occurred 

s i 

N1 

Field Intensity 

:ft) 

RDN or FLF 

Day Allowable Dose Occurred 

N 2 

N’2 

Probability of Fatality 

F l 

PI 

Probability of Noneffective 

P 2 

P2 


In the input sections of the program, if a programing symbol given in 
Table 3-1 is followed by L, H, or S, the ourrbvr represented by the symbol is 
either the lowest, highest, 01 Incremental valu taken by the variable repre¬ 
sented by the preceding letters TAL means the lowest value the time of 

arrival will rake), , , 








^* :«:rwrzt* are W/ed t-: •* <- /Uw diagram by the or|r.i ratter cf 
• tat raw t lire pva^et with the . Specifically, each of the 

coding lire numbers art four dig: " irst digit represents the major routine 

being calculated («ee Table B-II); second and third digits correspond to the 
flew diagram box number (xx) ; the fourth or unlts-posltlon digit Is "open” to 
allow for sequencing corrections, additions, and expansions (0). 

TABLE B-II 

. Major Program Routines ______ 

Code 

Number Series 


Input 

lxxO 



Calculation of Buildup 


Duration in hours 

3xx0 

Initial Conditions 

4xx0 

Calculate 1st Day Doses 

5xx0 

Day by Day Doses 

6xx0 

Outcome Probabilities 

7xx0 

Formats & ,f Housekeeping" 


Steps 

8xx0 


2. Explanation of the Flow Diagrams 

The first section of the program (block numbers 1-36) initializes the 
input data and sets up the control loops for the variables. Block numbers 36 
through 40 select the appropriate duration multipliers to calculate the dura¬ 
tion of Particle Buildup (T3). 

T3 - (1.133 + ETA)*TA when 2 =1. 

---—_ ( 1.175 + ETA)*TA-when 2 - 1.2-.- (B-13) 

■ (1.200 + ETA)*TA when 2 - 1.4 

the range of Buildup Duration Factor (ETA) is between zero (corresponding to 
the displaced sine function) and five (corresponding • ‘th 5TA ^ - TA function). 

The program sets the initial conditions in block 47. Clock Time (T) begins 
at time of arrival (TA) and the Day Counter (N) begins at one, the time is 
incremented by 3 minutes or 0.05 hrs. 
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Field ir.ttr*ity (RDN) ter unit of tine is calculated ir. MocVs 4i, 45. 
usin* the product of Reference Irtensltv (r/hr). Particle Arrival functlcr, and 
Decay function. 

Limit on T3 is: 

RDN - | Jl - cos(2^)«rjT* Z 0 < < - (B-14) 

The test, block 42, determines if there is any movement outside in fallout--or 
if T2 is positive. 

Block 52 is similar to 47 except the field is attenuated by the protection 
factor. Block 53 tests to determine if buildup of particles is continuing. 
Block 54 increments time out of buildup phase so that field intensity is 
governed only by the decay function corresponding to the algebraic formula. 

0 - T“ Z . (B-15) 


Block 56 tests to determine if the 1st day is completed. Block 55 increments 
time by 12 minutes during buildup. Blocks 49, 50, 51 are identical to 44, 45, 
46 but are not linked because of a change in the time intervals. 


FLD in block 59 is equivalent to RDN but is used to account for dose when 
the time step is by days. 


r(t) 


r 24 + TA t -Z dt 
T 


1-Z|25+TA 

C(t) - 1 ^ ■ | T when Z 4 1 (B-16) 

f(t) „ 


Block 60 does not have a time term since the integration in block 59 was over 
a finite time. 


Blocks 64, 65, 66, 67 and 68 sets up transition from the first day, to 
day-by-day calculations. At blocks o4-C3 " u '' t<rtc ERD division into recoverabl 
and nonrecoverable fractions is made. Block 66 steps time by 12 hours. Block 67 
sets maximum dose register. Block 69 solves equation for Field Intensity summed 
over 24 hours but corrected according to PF, 


FLD 


R»T~ Z *(24 hrs.) 
PF(K) 


(B-l7) 
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72 .pit let tit recovercblo dot*. Biock 71 th t nert ecoverai It 

dose. Block. 72 teats for sure death. Block 74 test* the day'* do*e vi. )Uxla.r. 
ERD. If there la a new Maxirun ERD, it 1* *et in block 75, and the day and tine 
•re incremented In blocks 76 and 77. If t K e peak ERD (D3) has been passed 
(D < D3), then the next point of interest is whether or not the allowable level 
has been reached. This test is block 78. If the dose is above the allowable, 
block 79 tests to see if the allowable is reachable--i.e., nonrecoverable dose 
is less than the allowable. If the nonrecoverable dose exceeds the allowable, 
then the program is in a loop as the length of stay in shelter would approach 
infinity—therefore, the error message. 

Blocks 80 through 92 calculate the probability of fatality and the probabili 
ty of noneffectiveness. 

The preceding explanation should give some understanding of the computer 
program and how it represents the mathematical formulation. This program, then, 
generates the possible data combinations that are analyzed in this volume. 








j 

i list ~ULLL*11 1 / 



Code Kumber Series 1XX0 
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PROGRAM MAINLINE 

01 MENS I ON PFNtlO) „ __ _ . .. 

N2*N2 

KO0E?0 ..____ ____ 

READ INPUT TAPE 5, IOOO.RL,RH,RS 

_ READ INPUT TAPE 5. 1000.TAL.TAH.TAS _ 

READ INPUT TAPE 5,1000,T2L,T2H,T2S 

READ.INPUT TAPE 5,1OOOj Q4L i D4H* f)4J>_ __ 

READ INPUT TAPE 5.1000,ZL.ZH,ZS 

READ_JNPUT TAPE 5. 1000. FL . FH,F S___ _ 

READ INPUT TAPE 5* 1000 * 8L 5 8H* 8S 

_ READ INPUT TAPE 5. 1 QUO. ETAL * ETAH. ETAS _;__ 

READ INPUT TAPE 5,1001,PFN111,PFNI 2 I *PFNI 3 I,PEN(4),PFN(51,PFNI 6), 

IPFNI n,PENm,PFNj9>__ 

READ INPUT TAPE 5,1002,KMAX 

1000 FORMAT ( 3F15 .5i)_1_1_ 

1001 FORMAT(9F6.0) 

1002 FORMAT! 110) ___ 

WRITE OUTPUT TAPE 6,8933 

8933 FORMAT ( 115H_R_TA_I2__04._Z __F 0 

1 ETA K Pi P 2 D3 Nl 0 N 02 1 

2110 _R=RLrWS __:..... 

2120 IF1R-RH12130,2130,8820 

2130 H*R + RS __ 

2140 TA = T AL-TAS 

2150 IFITA-TAH12l60j2160,2120___ _ 

2160 TA= T A + f AS 

2170 _T2 = T2L-T2S _..._____ 

2180 IFlT2-T2H)2190,2190,2150 

2190 T2 = T2»T2S ___ 

2200 D4-04L—D4S 

2210 IF(04-D4H)2220,2220,2180_ 

2220 D4*D4+04S 

2230 Z*ZL-ZS _ ____ 

2240 IF(Z-ZHI2250,2250,2210 

2250 Z*Z + Z S __ 

F=FL-FS 

22 70_ I F (F-FH) 228 0,2280,22 40__.___ 

2280 F=F+fS 

' ■_B a BL—BS_j___ 

2300 1F(B-BHJ2310,2310,2270 

2310 B-B+BS _ 

ETA=ETAL-FTAS 

2330 IFIETA-ETAH1234 0,234 0 ,2300 __ 

2340 ET A* ETA + ETAS 

2350 DO 9000 K*1,KMAX_____ 

3360 IFIZ-l.1)3370,3380,3380 

3370 T3*(1.1333+ETA)»TA __ 

GO TO 4000 

_3380 IF1Z-1.1313390,3400, 3400 __ 

3390 f 3 = (1•175+ETA1 •TA 

GO T 0 4000 ___ 

" 3400 T3*I1.20*ETA)*TA 

4000 CONTINUE _ 

4410 D*0• 

Dl-O. _ _____ 
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4420 

T-TA 

IFIT-IT2+TA)14430,5480,5480 


4430 

T-U0.05 


4440 

RDN-IT-TA)*3.i4l6/T3 


4450 

RON-l.O-COSFI RON) 


4460 

RON* CR-RDN)•IT--ZZ)/2.0 


4470 

0*0+IRDN-0.05) 



GO TO 4420 


5480 

T-T+O.l 


5490 

RON-(T-TA)*3«1416/T3 


5500 

RON-1.O-COSFIRDNI 


5510 

RDN*IR»RDN)'»(T-*ZZ>/2.0 


5520 

0*0*(R0N-0.2/PFNI K )) 


5530 

IF C T-ITA+T3)>5550,5550,5540 


5550 

T-T*0.2 

GO TO 5490 

*•' 

5540 

T-T+O.l 

IFIT-ITA+24.0)J5570.5610.5610 


5570 

IF(Z-1.015590,5580,5590 


5580 

ARG-24.0/T 



FLO-LOGFIARG) 

GO TO 5600 


5590 

ZZZ-l.O-Z 

FLO* I <T*»ZZZ)-UTA*24.0)**ZZZ))/<Z-1.0) 


5600 

D-O+IFLO-R/PFNIK)) 


5610 

IF 10-700.16000,9000,9000 


6000 

CONTINUE 


6640 

Dl-F-D 


6650 

02-11.0-FI*0 



66 70 
6680 


6750 

6760 


6780 

6790 


7710 

7000 


7800 

■781 Q 


7820 

7830 


D3-0 

N»N»1 


6690 

6700 


6710 

J02SL 


FL0-IR-24 .OlMT — ZZl/PFNIK) 
01»I l.-B)r01»IF»FLP) _ 


02-02 

P»P1 


2 * 1 1.0 
♦02_ L 


f-F) »FL0 


6730 

6740 


IF(0-700.16740*9000,9000 
IF ID-03>6 780.6750,6 750 
03-0 


Nl-N 


6770 


T«T*24.0 
GO TO 6680 


IF I 0-0417000,6790,6790 
IFID2-04>6770,7710,7710 


KODE-1 

CONTINUE 


N2-N 

1FIQ3-1QQ,>9000,7820,7820 


IF 103-200.) 7880*78 30,7830 
IF I 03-300,I 7890,7840,7840 














7840 IF(03-700.)7900,9000,9000 

7880 P2 = (03- 100 . ) /200 . _ _ ___ ______ 

7910 Pl=0. 

CO TO 8000 _ . . . ... 

7890 P2=((J.O*D3>-100.1/1000. 

_ CO TO 7920 __ 

7900 P2=(700 .-031/500. 

7920 Pl*(03-200. 1/500. _ ___ __ _, 

8000 CONTINUE 

6950 MR I T 6 OUTPUT TAPE 6,6934jKODE*R,TA,T2.04,2,F,B,ETA,PFN(K1, 
1PI,P2,D3,NI,0,N2,D2 

8934 FORMAT(I4.4X.F6.0.FB.1.F8.J.F6.0.F6.2.2X.F5.3.2X.F5.3.2X.F5 
1F5.0,F7.3,F7.3.F7.1,I5,F7.1,I5,F9.2> 

9000 CONTINUE _ _____ _ . 

K0DE*0 

__ _G0 I0_2330 __' 

8820 CONTINUE 

" '.END SENSPROG ~ .. " — —: ~ ; ~ - ;--:: 


END 
F INI S 


300. 

3000. 

1200. 

1.0_ 

6.0 

3.0 _ 

0.0 

.5 

.3 

20. 

70. 

30. 

1. 

1.3 

.2 

.85 

.94 

.05 

.02 

.029 

.005 

0 . 

5.0 

2.825 

10. 

40. 100. 

500. 
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